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ABSTRACT 



Context. Observations of quasars shining through foreground galaxies, offer a way to probe the dust extinction curves of distant 
galaxies. Interesting objects for this study are found in strong gravitational lensing systems, where the foreground galaxies generate 
multiple images. 

Aims. The reddening law of lensing galaxies is investigated by studying the colours of gravitationally-lensed quasars, and a handful 
of other quasars where a foreground galaxy is detected. 

Methods. We compare the observed colours of quasars reported in the literature, with spectral templates reddened by different extinc- 
tion laws and dust properties. The data consists of 21 quasar-galaxy systems, with a total of 48 images. The galaxies, which are both 
early- and late-type, have redshifts in the interval z = 0.04 - 1.51. 

Results. We measure a difference in rest-frame B - V between the quasar images we study, and quasars without resolved foreground 
galaxies. This difference in colour is indicative of significant dust extinction in the intervening galaxy. Good fits to standard extinction 
laws were found for 22 of the images, corresponding to 13 different galaxies. Our fits imply a wide range of possible values for 
the total-to-selective extinction ratio, Ry. The distribution was found to be broad with a weighted mode of ^v=2.4 and a FWHM of 
ARv=2.7 {(TKy ~ 1.1). Thus the bulk of the galaxies for which good reddening fits could be derived, have dust properties compatible 
with the Milky Way value (Ry = 3.1). 

Key words. dust,extinction - galaxies: ISM - galaxies: evolution - quasars: general - techniques: photometric 



1. Introduction 

Measuring the reddening due to dust in galaxies, provides infor- 
mation on dust grains in the interstellar medium. It is in addition 
important for dimming corrections, when measuring cosmologi- 
cal distances using Type la supernovae, the most direct probe of 
the expansion history of the Universe. 

It is often assumed that the wavelength dependence of ex- 
tinction in distant galaxies is identical to the average in the Milky 
Way. Ho wever, there are indications that this is not always the 
case (e.g. iPrevot et al.lll984t ICalzetti et alj|2000l) . The purpose 
of this work is to examine the reddening relations in galaxies 
over a wide range of redshifts, to improve our understanding of 
galaxy dust properties. 

For galaxies in which individual stars are resolved, dust ex- 
tinction can be studied using the pair method. In this method, 
the galaxy extinction curve is determined by comparing the 
spectral energy distribution of reddened and non-re ddened stars 
of th e same spectral type and luminosity class (Mas sa et alj 
I1983-) . For more distant galaxies, other methods must be used, 
such as comparing the galaxy colours with a dust-free model 
(although this is subject to large model dependencies) (see 
e.g. Patil et al. 2007), determini ng the differential extinction of 
multiply-imaged quasars (e.g. iFalco et aTl Il999l) . or measur- 
in^ the absoliite ext inction by foreground galaxies of quasars 
JOstman et alj|2006l) . There are problems using the method of 
differential extinction for multiply-imaged quasars. To obtain a 
unique Ry value, one of the images must have negligible dust 
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extinction when compared to the other, or the extinction curves 
must be similar along the two separate lines of si ght. The po- 
tential probl em of this hmitation is discussed in Mc Gough et al.l 
(l2005h and lEhasdottir et al.l (l2006l) . If these conditions are not 
met, the Ry value derived will be a function of the extinctions 
al ong the two lines of s ight, and their Ry values (see Eq. 6 
in 'McGough et al.l (l2005h ). In this paper, this restriction is cir- 
cumvented by studying the individual colours of background 
quasar images in order to investigate the reddening by fore- 
ground galaxies. This is possible because of the fairly homoge- 
neous spectral pro perties of qua s ars. T he method is very similar 
to the one used in Ost man et al.l ([2006), where the focus was on 
finding systems by matching the positions of quasars with galax- 
ies. In this paper, the focus is on using systems found through 
gravitational lensing. 

Studies of colours of Type la supernovae have also been 
used to constrain the dust properties of their host galaxies 
(e.g. 'JSeveer 1983; Branch & Tammann 1992: Riess et al.lll996l; 
iReindlet al...2005: .Guv et al... 20051 When the reddening cor- 
rection has been optimised to minimise the scatter of the Type la 
supernova Hubble diagram, a quite different wavelength depen- 
dence of the extinction curve, when comp ared to the Milky Way, 
has been derived (e.g. lAstier et aLll2006h . It is unclear whether 
these findings indicate the presence of non-standard dust along 
the line of sight to supernovae, or if they point to a different 
source of the brightness-colour relation, perhaps intrinsic to su- 
pernovae themselves. We aim to establish if the reddening laws 
derived using supernovae are applicable to different astrophysi- 
cal scenarios, such as quasars shining through galaxies. 
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In Sec.|2] different parameterisations of dust attenuation are 
discussed. The method for estimating the dust properties is pre- 
sented in Sec. [3] The resuhs and a discussion of the results are 
given in Sec.|4] Finally, we summarise and conclude in Sec.|5] 



2. Dust extinction parameterisations 

Extinction curves are often parameterised by the reddening pa- 
rameter called the total-to-selective extinction ratio Ry, which is 
defined by 



Rv = 



EiB-V) 



(1) 



where Ay is the V-band extinction and E{B - V) is the colour 

excess. 



E(B -V) = (B- y)obs -(B- V)u 



(2) 



The observed colour is denoted by (B - y)obs and the intrin- 
sic colour by {B - y)intr- The reddening parameter is a mea- 
sure of the slope of the extinction curve, and is related to the 
size of the dust grains. Large grains produce greyer extinction, 
with larger values of Ry- The average value of Ry in the Milky 
Way is 3.1, but the value varies wi thin the inter val 2 < Ry < 6 
between different regions (see e.g.'Drainel l2003L and references 
therein). The variation s are rare an d appe ar to be conn ected with 
dense molecular gas jjenniskens & Greenberd[l993h . Among 
distant galaxies, the vari ation of Ry could be larger. For ex- 
ample, iFalco et al] d 19991) derived values in the interval 1 .5 < 
Ry < 7.2 using the method of differential extinction between 
images of gravitationally-lensed quasars. Measurements of the 
mean Ry from Type la supernova samples yield lower values 

than the Milky Way a verage. Some examp les are Ry ~ 1-8 

( lKrisciunaset al.l l2000h. Rv = 2.5 5 + 0.30 (iRiess etal.lll"996h. 
Ry ^ 2.6 + 0.4 d^ilHps et al.' 1999^, Ry = 1.09 dTrippI 
119981). Ry = 2.5 (JAl tavilla et al. 2004) , Ry = 1.75 + 0.27 
JNobili & Goo bar 200 7J) and P = 1 . 77 ± 0.16 corre sponding to 



Ry = 0.77 + 0.16 ( Astier et al.ll2006l:lGuv et £1112007;) . Although 
parameterised as an interstellar extinction law, it is possible that 
the brightness-colour relation observed in Type la supernovae, 
could be intrinsic to the supernovae or related to circumstellar 
dust. 

The extinction curve of gala xies is ofte n described using the 
parameterisations by either Car delh et al.l (1989) (CCM), or by 
iFitzpatrickl (1 19991) . However, some galaxies, such as the Small 
Magellanic Cloud (SMC), are not well described by these laws. 
An alternative extincti on law, valid for the SMC, was presented 
bv iPrevot et al.l (Il984l) . Starburst galaxies are believed to have 
a different ex t inction law again, which has been described by 
ICalzettietal] ( |2000|) (CAB). All of these extinction laws are 
shown in Figure [T] 

Since the extinction curves for the Small Magellanic Cloud 
and starburst galaxies look different compared to the average one 
for the Milky Way, the extinction curve of the Milky Way may 
not be representative of high-redshift galaxies. Yet, it is often 
assumed that the Galactic laws can be applied to dust in other 
galaxies, e.g. when correcting the light from supernovae for host 
galaxy extinction. In this paper, such an assumption is relaxed. 



3. Fitting dust extinction properties 

We have gathered data for a small number of quasars with fore- 
ground galaxies, from the literature. We require that the red- 
shift of both the quasar and the galaxy have been determined 




1/A i^Lm-') 

Fig. 1. The CCM parameterisation of the Milky Way extinction 
curve (solid black line) for four different Ry (0.5, 2, 3.1 and 4) 
where the highest values of Ry correspond to the flattest curves. 
We also show the Fitzpatrick parameterisation (dotted blue) with 
Ry - 3.1, the Prevot law for SMC-like extinction (dashed green) 
with Ry = 3.1 and the starburst extinction law (dashed dotted 
yellow) with the pref erred value for starburst galaxies Ry - 4.05 
(ICalzetti et al.ll2000h . A^ is the extinction at wavelength A and 
Ay is the visual extinction. 



spectroscopically. We only consider systems for which appar- 
ent magnitudes have been measured in at least four different fil- 
ters, which are in the wavelength range of the redshifted, quasar 
spectral template. Given these restrictions, our sample of quasar- 
galaxy systems consists of 21 foreground galaxies with a total 
of 48 quasar images. Information about these pairs are listed in 
Table \T\ Most of these are found from strong lensing surveys. 
Ho wever, three of the m comes from analysing the quasar spec- 
tra (IWang etaTI 12004) and an additional three are found from 
coordinate matching JOstman et al.ll2006l) . 

To test whether our sample of 21 quasars with foreground 
galaxies differs in the colours from the main bulk of quasars, 
we estimate the rest-frame B - V colour, by K-correcting from 
the observed bands, for our set of quasars and compare it with a 
referenc e set from release three of the Sloan Digital Sky Survey 
(SDSS) (i Schneider et al.ll2005h . Table |2] contains the estimated 
values of B-V for our images, and Figure|2]shows the histogram 
of the two distributions. We find that our set is redder than the 
comparison set. The difference in mean rest-frame B-V be- 
tween the two sets is 0.7 magnitudes. There is also one quasar, 
WFI2033-4723, for which three out of its four images are signifi- 
cantly bluer than the reference set. The two distributions of rest- 
frame colours were compared using the Kolmogorov-Smirnov 
test, yielding a negligible probability that the difference is only 
due to statistical fluctuations. We conclude that there must be 
quasars in our sample whose colours are affected by reddening 
in the intervening galaxy. 

The quasar spectral template use d, is a combinatio n of the 
HST radio-quiet composite spectrum (Telfer et al. 20 02|), and th e 
SDSS median composite spectrum ( Vand en Berk et al.l 1200 lb . 
To simulate the effects of dust attenuation, the template spec- 
trum is redshifted and reddened, by different values of Ry and 
E{B - V). This enables synthetic colours to be estimated for dust- 
reddened quasars. To determine the values of the dust parameters 
that best describe the observed magnitudes, we compare the ob- 
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Table 1. Information about the quasar-galaxy pairs in our sample, zq is the redshift of the quasar and zc is the redshift of the galaxy. 
The last two columns list what filters were used in our analysis and the references for the magnitudes. 



Quasar 


Images 


Zq 


Zc, 


Galaxy type 


Filters 


Source" 


SDSSJ131058.13+010822.2 


A 


1.39 


0.04 


star forming 


ugriz 


1 


Q2237+030 


ABCD 


1.69 


0.04 


late 


F160W F205W F555W F675W F814W 


2 


SDSSJ114719.89+522923.1 


A 


1.99 


0.05 


starburst 


ugriz 


1 


SDSS J084957.97+5 10829.0 


A 


0.58 


0.07 


starburst 


ugriz 


1 


SDSS 1155+6346 


AB 


2.89 


0.18 


early 


F160WF555WF814WK 


2,3 


MG 1654+ 1346 


A 


1.74 


0.25 


elliptical 


F160W F555W F675W F814W 


2 


CXOCY J220132.8-320144 


AB 


3.90 


0.32 


spiral 


griKs 


4 


SDSS J0903+5028 


AB 


3.61 


0.39 


early 


F160WF555WF814Wri 


2,5 


SDSS J0924+0219 


ABC 


1.52 


0.39 


elliptical 


u g r i F160W F555W F814W 


2,6 


CLASS B1152+199 


AB 


1.02 


0.44 


late 


F160W F555W F814W I R V 


2,7 


HE0435-1223 


ABCD 


1.69 


0.45 


SO 


F160WF555WF814Wigr 


2,8 


Q0142-100 


AB 


2.72 


0.49 


early 


F160W F555W F675W F814W 


2 


HE0230-2130 


Al A2 B C 


2.16 


0.52 


SO/Sa 


BRIKF814WF555W 


2,9 


BRI0952-0115 


AB 


4.43 


0.63 


elliptical 


F160W F555W F675W F814W 


2 


WFI2033-4723 


Al A2 B C 


1.66 


0.66 


Sb/Sc 


F160WF555WF814Wi 


2, 10 


SDSSJ1004+4112 


ABCD 


1.74 


0.68 


cluster of galaxies u g r i z F160W F555W F814W 


2, 11 


J 1004+ 1229 


AB 


2.65 


0.95 


elliptical 


F160WF555WF814WI 


2 


MG0414+0534 


Al A2 B C 


2.64 


0.96 


early 


F160W Fl low F205W F675W F814W 


2 


SDSS J012147.73+002718.7 


A 


2.22 


1.39 


unknown 


ugriz 


1 


SDSS J145907.19+002401.2 


A 


3.01 


1.39 


unknown 


ugriz 


1 


SDSSJ144612.98+035154.4 


A 


1.95 


1.51 


unknown 


ugriz 


1 




° 1 : ' Adelman-McCarthv & for thf 


; SDSS Collaboration 
.(:2003), 7::Toftetal, 


('2007), 
200Q), 


2:'Kochanek et a 
8:iWisotzkietal. ( 


.' ('2007), 3:'Pindor et al." ('2004"), 4:'Castander et al.' ('2006) 


Johnston et al. (2003), 6: Jnada et al 


2002), 9: Wisotzki et al, (, 1999,), 10: .Morgan et al., (,2004,), 


Oauri et al.l d2004l) 
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Table 2. Estimated rest-frame B 
quasar images in our sample. 



V colour for the different 



Quasar 


Images 


B-V 


SDSSJ131058.13+010822.2 


A 


0.8 


Q2237+030 


ABCD 


3.4,2.1,3.6,3.4 


SDSS J084957.97+5 10829.0 


A 


1.1 


SDSS 1155+6346 


AB 


1.5, 1.1 


MG 1654+ 1346 


A 


0.4 


CXOCY J220132.8-320144 


AB 


1.7, 1.8 


SDSS J0903+5028 


AB 


0.7, 0.5 


SDSS J0924+0219 


ABC 


0.4, 0.6, 0.2 


CLASS Bl 152+199 


AB 


1.8, 1.7 


HE0435-1223 


ABCD 


1.1, 1.0, 1.0, 1.0 


Q0142-100 


AB 


0.1,0.1 


HE0230-2130 


Al A2 B C 


0.4, 0.2, 0.3, 0.8 


BRI0952-0115 


AB 


0.7, 0.8 


WFI2033-4723 


Al A2 B C 


0.4,-0.1,-0.3,-0.1 


SDSSJ1004+4112 


ABCD 


0.1,0.6,0.3, 1.6 


J 1004+ 1229 


AB 


3.9, 5.0 


MG04 14+05 34 


Al A2 B C 


2.9,3.1,2.8,2.8 


SDSS J012147.73+002718.7 


A 


1.0 


SDSS J145907.19+002401.2 


A 


0.9 


SDSS J144612.98+035 154.4 


A 


0.8 



served colours with the synthetic colours for difi'erent Ry and 
E{B - V). When choosing colours X - i - j to compare among 
all possible filter combinations, we want to use the colour com- 
binations that have the greatest potential to separate cases of red- 
dening from those of no reddening. We select the combinations 
that provide a large difference between observed colour and ex- 
pected colour without dust extinction, E{i - j), and have small 



cr 2> 




Fig. 2. A histogram over rest-frame B-V for the quasars in 
SDSS DR3 (solid line), and for the quasar images in our sam- 
ple (dashed line). 



errors. We select colours by maximising the sum over Qij, which 
is defined by 



Qi 



[Eii-j)f 
4 



(3) 



where cTij is the total error of the colour difference, i.e. the er- 
rors of the observed magnitudes combined with the uncertainty 
in the template colours. The sum over Qij runs over /- 1 colours 
where / is the number of filters, and each filter must be used at 
least once. The colours chosen from maximising the sum over 
Qij are used when calculating the x^ described below. For exam- 
ple, when observations are available in the ugri filter set, the six 
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Fig. 3. Uncertainties in the template colours for u-g,u-r,u-i, 
u - z, g - r, g - i, g - z, r - i, r - z, i - z (from black to light grey) 
as a function of redshift z- 

possible colour combinations are u - g, u - r, u - i, g - r, g - i 
and r - i. We assume that Qug and Qu, have the largest values of 
Qij. From the remaining four combinations, u-r, g-r, g-i and 
r - /, we can only choose from u- r, g- r, and r - /, because the 
colour M - / can be constructed by combining the colours we have 
already picked, u- g and u - i. We then choose the filter combi- 
nation, of the three remaining, with the largest Qij to obtain the 
three colours to be used in the analysis. 

For each quasar image, the best-fit dust extinction is 
determined by comparing the synthetic colours, Xsyn - 
^syn[^v, E(B-V)], with the observed ones, Xobs, and minimising 
the;if2 = X^[Rv, E(B - V)], which is defined to be 



X^ = (^s 



syn 



■Xobs)^^"'(Xs 



,yn 



■ ^obs)- 



(4) 



"V is the covariance matrix including the errors from the ob- 
served and synthetic colours. The errors of the synthetic colours 
account for both the intrinsic variability of quasars, and the pos- 
sible dust extinction in the host galaxy of the quasar. These er- 
rors are calculated from a set of quasars with similar redshifts 
(|Az| < 0.05), without any resolved foreground galaxies. These 
were taken from the third data release of the Sloan Digital Sky 
Survey ( Schn eider et al . 2005). Their colours were /T-corrected 
from the SDSS filters, to the filters used in the observations con- 
sidered in this work. To avoid a large impact from colour out- 
liers, quasars with a colour differing by more than three stan- 
dard deviations from the mean were excluded. About 12% of the 
quasars were rejected by this cut. In Figure |3] the uncertainties 
for the colours from the template are shown for the SDSS filters 
at different redshifts. As can be seen in the figure, all colours 
involving the u band have very large errors for quasars with 
z > 2.5. It is also at z ~ 2.5 that the number of available compari- 
son quasars falls drastically. The covariance matrix also includes 
the correlations between estimated magnitudes in diff'erent filters 
and consideration to the fact that the same filter can appear more 
than once in a set of colour comparisons. 

The quasar images were fitted with both the CCM law and 
the Fitzpatrick law with Ry and E(B - V) as free parameters. 
We also considered SMC like extinction, according to the Prevot 
law. The CCM law that we use is the improved version of the 
relation by Cardelli et al. (1989), where the extinction law, for 
the wavelength range 1.1 yum ' < 1/A < 3.3 um~', has been ex- 
changed for the curve provided bv IO' Donnelll i 1 9941) . The Prevot 



law used here consists of tw o parts. For A < 27 00 A, we use a 
linear fit to the SMC data in lPrevot efal] (!l984'). and at longer 
wavelengths we use the extinction law by|Fitzpatrick (1999). 

In addition to dust in the foreground galaxy, we have also 
investigated the possibility of extinction in the host galaxy of the 
quasar. 



4. Results and discussion 

Of the 48 images, 26 yielded inconclusive results. The chi- 
square of their fits was so high that, if they had been affected 
by dust following our extinction laws, there was only a 0.5% 
probability that such a high chi-square value would be obtained. 
However, a poor fit does not automatically imply that there is 
a strange, unknown reddening law. It can also be the result of 
for example (1) underestimated observational errors, (2) the ex- 
tinction law having more degrees of freedom than the laws we 
have assumed, (3) peculiar quasars that have an intrinsic spec- 
trum that differs significantly from the template, (4) effects from 
other intervening objects that are unresolved in the available im- 
ages, (5) microlensing effects, or (6) changes in colour due to 
the time variability of the quasar. The significance of the latter 
t wo possibilities, when c omparing quasar images, are discussed 
in lYonehara et all (l2007l ). 

The fit results for the remaining 22 images (corresponding 
to 13 galaxies) are presented in Table [3] In Figure |4] a few of 
the confidence levels obtained from the x^ fitting of the quasar 
images are shown. The error intervals in Table |3]indicate the Icr 
uncertainty. The errors are calculated with the assumption that 
the extinction law is valid for a positive colour excess, and for 
all positive values of Ry up to 6.0. Negative values of E(B - V) 
are not expected for normal quasars, affected by dust. However, 
an intrinsically blue quasar is best fitted with a negative E{B-V). 
For small and negative values of Ry (< 0.7), the extinction laws 
have a peculiar behaviour (see e.g. the CCM parameterisation 
with Ry = 0.5 in Figure [T]i. The non-smooth function, and the 
negative values of the extinction A^/Ay, are probably hard to ex- 
plain using a physical dust model. We chose however to consider 
all fits, down to Ry - 0. For large (positive) values of Ry, the 
extinction becomes less wavelength-dependent, and it becomes 
more difficult to distinguish from the case of no extinction. When 
the colours of a quasar image in our analysis can be explained 
without dust, there is no limitation on Ry, because there is no Ry 
dependence in the extinction laws for E(B - V) = 0. 

The law which we refer to as the Prevot law, consists of two 
parts. For A < 2700 A, we use a linear fit to the SMC data in 
iPrevot et al] ([19840, and for A > 2700 A, we use the parame- 
terisation of the Milky Way extinction by Fitzpatrick ( 1999). If 
the foreground galaxy, where the extinction occurs, is located at 
low redshift, most of the filters used for observations will corre- 
spond to regions of the spectrum that are affected by the second 
part of the extinction curve, which is identical to the Fitzpatrick 
law {A > 2700 A). If, on the other hand, the extinction redshift 
is high, the filters will correspond to regions of the spectrum that 
are affected by the first part of the extinction curve. This part has 
the relation A,i/£'(B-y) = a -i-/7x-i-/?y. It is therefore not possible 
to estimate Ry from the observed colour, e.g. B -V, because the 
Ry dependence is subtracted out, B-V - As-Ay = E{B-V){a + 
bxB + Ry) - E{B - V){a + bxy + Ry) = E(B - V){bxB - bxy). This 
is not the case for e.g. the Fitzpatrick or CCM law. For objects 
which have high extinction redshifts, no preferred value of Ry 
can be determined for SMC-like extinction, using this parame- 
terisation. The same problem occurs for the CAB law, where a 
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Table 3. Best fits of Ry and E{B - V) together with the 1 cr uncertainties. Note that the lowest possible value for Ry here is 0. 1 , since 
Ry = usually is not defined. Also note that for some redshifts and filter combinations, there is no Ry dependence for the Prevot 
law. The interval is then given as 0.1 - 6.0. Zd is the redshift of the dust extinction, which is zc for foreground extinction and zq 
for host extinction, d is the impact parameter expressed in kpc. When this was unavailable, it is marked with a hyphen in the table. 
When several laws and/or dust redshifts is either fitted for one image, they are all listed in the table below. The probability for each 
to obtain a higher ;^f^ value than what was obtained for the fit is given in the last column. We include all fits with a probability larger 
than 0.5%. The different images are separated by a horizontal line. 



Quasar 


Image 


Zd 


d (kpc) 


Rv 


E(B - V) 


Law 


Prob 


SDSSJ131058.13+010822.2 


A 


0.04 


15.0 


3.5 (1.4-6.0) 


0.3 (0.2-0.4) 


CCM 


96% 


SDSSJ131058.13+010822.2 


A 


0.04 


15.0 


3.5 (1.8-6.0) 


0.3 (0.2-0.4) 


Fitzpatrick 


94% 


SDSSJ131058.13+010822.2 


A 


0.04 


15.0 


3.5 (1.8-6.0) 


0.3 (0.2-0.4) 


SMC 


95% 


SDSSJ131058.13+010822.2 


A 


1.39 


15.0 


0.1 (0.1-2.5) 


0.1 (0.1-0.2) 


Host CCM 


12% 


SDSSJ131058.13+010822.2 


A 


1.39 


15.0 


1.1 (0.3-2.2) 


0.1(0.0-0.2) 


Host Fitzpatrick 


46% 


SDSSJ131058.13+010822.2 


A 


1.39 


15.0 


0.1 (0.1-6.0) 


0.1(0.1-0.2) 


Host SMC 


20% 


Q2237+030 


A 


0.04 


0.7 


0.1 (0.1-0.2) 


0.7 (0.5-0.9) 


Fitzpatrick 


1 % 


Q2237+030 


A 


0.04 


0.7 


0.1 (0.1-0.2) 


0.7 (0.5-0.8) 


SMC 


2% 


Q2237+030 


C 


0.04 


0.6 


0.1 (0.1-0.2) 


1.2(0.9-1.4) 


CCM 


11 % 


Q2237+030 


C 


0.04 


0.6 


0.2(0.1-0.6) 


0.8(0.6-1.0) 


Fitzpatrick 
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Fig. 4. Confidence levels for some of the quasar images that could be fitted by dust extinction in the intervening galaxy using the 
Fitzpatrick extinction law. The levels correspond to Icr for one parameter (black line) and 68% (dark blue region) and 90% (pale 
blue) for two parameters. 



preferred value of Ry cannot be determined for any extinction 
redshift, because the extinction law provides colours that are in- 
dependent of -Rv. This is because the parameterisation we use is 
AaIE{B-V) - a+bx+cx^+dx^+Rv, and thus B-V - A^-Ai/is 
not dependent on Ry- This should be remembered when studying 
the fitted values in Table [3] 

We note that even though we have considered both interven- 
ing and host extinction, only host extinction cannot explain the 
difference in B - V that we see, when comparing the colours 
of our quasars with quasars without intervening galaxies (see 



Figure |2]i. Furthermore, we have taken the possibility of dust ex- 
tinction within the host into account, in the calculation of the 
covariance matrix. 

In Figure |5] observed colours are plotted together with the 
best synthetic dust colours for each of the quasar images, that 
could be fitted by dust (i.e. the images in Table O. In the fig- 
ure, comparison is made with the synthetic magnitude at the 
longest wavelength arbitrarily adjusted to match the observa- 
tions. However, this was not the case in the actual fitting pro- 
cedure, which involves only differences in magnitudes. In all of 
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Fig. 5. Comparison of the measured magnitudes (black signs) with the fitted magnitudes (hght blue solid line) for the images that 
could be fitted using the dust laws. The dark blue dotted line shows the expected magnitudes without dust. Note that sometimes only 
one line is visible when the best fit is with little or no extinction. The ground-based measurements are marked with plus signs while 
the HST observations are marked with triangles. All magnitudes are normalised so that the magnitude at the longest wavelength is 
zero. 
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these cases, the fitted magnitudes agree very well with observed 
values. 



4.1. The badly fitted quasar images 

In Figure |6] we show how well the best extinction law agrees 
with observed magnitudes for each of the quasar images that 
could not be fitted with regular dust extinction. 

A possible explanation for the difficulty in fitting some of 
the quasar images with the reddened template is that in some 
cases, observations of the same image in different filters were 
made at different times. If the object is variable, and has changed 
significantly between the times of observations for the different 
filter sets (e.g. HST and ground-based), this could lead to erro- 
neous observed colours, which cannot be correctly fitted by dust. 
One possible example is SDSS J1004+4112 where the ground- 
based and HST observations appear to be shifted with respect 
to each other (see Figure |6l). Therefore, in the cases where we 
have enough data points, we separate the HST and ground-based 
observations to see if this will provide a good fit. There are 
1 1 quasar images that have observations in a minimum of four 
different filters, which are either ground-based or space-based. 
When applying our analysis to these data, we find that two of 
the images of HE0230-2130can be described by dust extinction. 
The best-fit values are provided in Tabled 

When calculating the errors of the synthetic colours, we de- 
fined colour outliers as quasars with a colour differing by more 
than three standard deviations from the mean. This led to the 
exclusion of 12% of the quasars in our reference set. Thus we 
would expect about 2.5 of our 21 quasars to be colour outliers, 
if our sample is as homogeneous as the SDSS sample. Some of 
the badly-fitted quasar data could correspond to so-called pecu- 
liar quasars, differing significantly from the template in one or 
several filters. 

If a deviation is observed in only one filter, it could be ex- 
plained by a strong absorption or emission feature, which is red- 
shifted into the wavelength region of the filter of observation. 
Some badly-fitted images, such as SDSS J144612.98H-035154.4, 
have colours that agree well with synthetic dust colours, apart 
from one band (see Figure |6]l. For these, we make the assump- 
tion that the diverging magnitude is erroneous and redo the fit 
without. Because observations in at least four filters are required, 
there is only one case of interest, SDSS J144612.98-I-035154.4. 
The best-fit values for this system, when removing the u band 
magnitude are provided in Table |4] 

4.2. Global results 

To understand how the fitted values of Ry are distributed, we 
assume that all quasars are affected by dust extinction in the in- 
tervening galaxy, following the Fitzpatrick law. For each image, 
we calculate a probability distribution for Ry, by locating the 
maximum ;^f^ probability for each Ry, when E(B - V) is set free. 
The maximum of the distribution will have the same value, as 
the probability for the Fitzpatrick law, in Table |3] The probabil- 
ity distributions for all images that could be fitted with dust, are 
then added together. In this way, we take into account the pos- 
sibility of non-Gaussian distributions, and the result is weighted 
by the probability of the fit. We exclude images for which our fits 
are consistent with almost no extinction because the Ry depen- 
dence is weak then. The result is presented in Figure |2l where it 
is evident that the inclusion of the objects in Table IH does not 
affect the shape of the curve significantly. Using the larger set. 




Fig. 7. The co-added probability distributions for the fits with in- 
tervening dust following the Fitzpatrick extinction law, weighted 
by the probability of the fit. The dotted line shows the added 
probability when using the well fitted quasar images in Table [3] 
and the solid line when also the objects in Table |4] are included. 
The added probabilities for the subset of early-type galaxies 
(dashed line), late-type galaxies (dash-dot) and starburst/star- 
forming galaxies (dash-dot-dot-dot) are also included in the fig- 
ure. 



we found a most probable value of Ry of 2.4, with a FWHM of 
2.7. This value is slightly lower than the Galactic value of 3.1, 
but well within the range due to the large spread of fitted values. 
In Figure[8] the best-fit values of Ry are plotted as a function of 
redshift and E(B - V). There are no clear trends in the depen- 
dence with these parameters on Ry. 

To test the possibility of a dependence of galaxy type for 
Ry, we divide Figure|7]into three plots, one for early-type galax- 
ies, one for late-type galaxies, and one for starburst/star-forming 
galaxies. These are shown in Figure|2] The plot for the late-type 
galaxies have less features than the one for all galaxies, but the 
most probable Ry remains similar. The plots for the early-type 
galaxies and the starburst galaxies are more different. However, 
these are constructed using only four and three images, respec- 
tively. It would be interesting to repeat the same study for a larger 
number of galaxies in the future. 

We note that since a large part of the galaxies studied are 
lensing galaxies, the sample will not necessarily be representa- 
tive of the universal galaxy population, because strong-lensing 
surveys are more efficient in finding massive galaxies. 



4.3. Specific systems 

Several of the quasar-galaxy systems that we study, have multi- 
ple images. Because the multiple images have the same intrin- 
sic colours, the difference in colour between the images should 
be explained by dust extinction, or possibly by microlensing or 
time-variability effects. In this section, we compare fits for dif- 
ferent images of the same quasar. We also compare with other 
dust property estima tes of the same object from the literature. In 
lOstman et alj ( l2006h . some of the objects here were investigated 
with a similar method, that has now bee n improv ed. Thejnost 
important improvement is that, while in lOstman et al.l (l2006l) . 
we compared all magnitudes with the value in the u band in the 
expression for the ;^f^, we here optimise the choice of colours 
to maximise the probability of detecting a potential dust extinc- 
tion (see Section|3j. Another improvement is that we include the 
observational errors, in addition to the template errors, in the co- 
variance matrix. Some of the systems that we have studied have 



L. Ostman et al.: Extinction properties of lensing galaxies 



Q2257 + 030 B, CCM, R.= 2.0 



SDSS1 155 + 6346 A, CCM, R,= 0.1 



SDSS J0924 + 0219 A, CCM, R,= 0.1 SDSS J0924 + 0219 B, CCM, R.= 0.1 




HE0435-1223 B, FTZ, R,= 0.1 




1/A [^m-) 
HE0230-2130 A2, CCM, R,= 6.0 




0.0 0.5 1.0 1.5 2.0 2.5 
1/i (,im-) 



WFI2033-4723 A2, CCM, R,= 0.1 



2.5 




- 


2.0 




y 


1.5 




A 


1.0 


/^ 


- 


0.5 




- 


0.0 




: 


5 







0.5 1.0 1.5 2.0 

1/X (Mm-) 

SDSS J1004 + 41 12 B, CCM, R,= 0.1 




0.5 1.0 1.5 2.0 2.5 5.0 
1/A (,2m-) 

J1004+1229 B, FTZ, R,= 0.8 

_ . A' \ 





0.0 0.5 1.0 1.5 2.0 

1/A (^m-) 

CLASS Bl 152+199 A, FTZ, R,= 0.1 



2.0 




- 


1.5 




/ 


1.0 




/A 


0.5 


y. 


+ 


0.0 


^ + 


- 


0.5 




- 



l/A (^m-) 

HE0435-1223 C, FTZ, R,=0.1 




HE0230-2130 B, CCM, R,= 6,0 

5 r 




0.0 0.5 1.0 1.5 2.0 2,5 
l/A (,im-) 



WFI2033-4723 B, CCM, R,= 0.1 




l/A (^m-) 

SDSS J1004 + 4112 C, CCM, R,= 0.1 




2.5 


^.^^ 


- 


2.0 


/^ N 


V 


1.5 


/^A+ 


\ 


1.0 

0.5 


A + 


+ - 


0.0 


A 


- 


n 5 







0.5 1.0 1.5 2.0 2.5 3.0 
1/X (,im-) 

CLASS Bl 152+199 B, FTZ, R,=0.1 




1.0 1.5 

lA (Mm-) 



HE0435-1223 D, FTZ, R,= 0.1 




0.5 1.0 1.5 2.0 2.5 

1/X (Mm-) 

HE0230-2130 C, CCM, R,= 4.6 




0.0 0.5 1.0 1.5 2.0 2.5 
1/X {^l^-) 



WFI2033-4723 C, CCM, R,= 0.1 




1/X (Mm-) 

SDSS J10Q4 + 4112 D, FTZ, R,=0.4 



+ 



1.5 2.0 2.5 3.0 
1/X (Mm-) 



SDSS J144612. 98 + 035154. 4 A, FTZ, R,= 1.2 



1.0 
0.5 
0.0 

n 5 


K 


A- 


+ 




0.5 


1.0 1.5 2.0 
1/X (Mm-) 


2.5 3.0 



+ 




1/X (Mm-) 



1.0 1,5 2.0 2.5 3.0 

1/X (Mm-) 



2.5 


^,^^ 


- 


2.0 


y^^ 


s. 


1.5 
1.0 
0.5 


/ A+ 

A + 


+ - 


0.0 


A 


- 


5 







0.5 1.0 1.5 2.0 2.5 3.0 
1/X (Mm-) 

HE0435-1223 A, FTZ, R,= 0.1 




HE0230-2130 Al , CCM, R,= 6.0 




3.0 0.5 1.0 1.5 
1/X (Mm-) 



WFI2033-4723 Al, FTZ, R,= 0.1 



^A 




0.5 1.0 1.5 2.0 

1/X (Mm-) 

SDSS J1004 + 4112 A, FTZ, R,= 0.2 



2.5 




^^ 


- 


2.0 




A ■ 


^^ 


1.5 




A 


•. 


1.0 
0.5 
0.0 


A 


+ + 


+ - 


-0 5 









0.5 1.0 1.5 2.0 2.5 3.0 

1/X (Mm-) 



J1004+1229 A, FTZ, R,= 1.2 




1/X (Mm-) 



Fig. 6. Comparison of the measured magnitudes (black signs) with the fitted magnitudes (light blue solid line) for the images 
that were badly fitted using the dust laws. The dark blue dotted line shows the expected magnitudes for no extinction. Note that 
sometimes only one line is visible when the best fit is with little or no extinction. The ground-based measurements are marked with 
plus signs while the HST observations are marked with triangles. All magnitudes are normalised so that the magnitude at the longest 
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Table 4. Additional best fits of Ry and E{B - V) together with the Icr uncertainties for quasar images which are likely to be dust 
extincted. The first three quasar images are badly described when using both HST and ground-based images, but are possible to fit 
with dust when using only the ground-based measurements. The last image is one where the u band magnitude seems to be deviant 
from the other magnitudes and we have made a fit without it. 



Quasar 



Image 



Zd 



d (kpc) 



Rv 



E(B - V) Law 



Prob Comment 



HE0230- 
HE0230- 
HE0230- 
HE0230- 
HE0230- 
HE0230- 



2130 
2130 
2130 
2130 
2130 
2130 



0.52 


6.1 


3.2(1.5-5.4) 


0.2(0.1-0.2) 


CCM 


13% 


ground 


0.52 


6.1 


3.4(1.8-6.0) 


0.2(0.1-0.2) 


Fitzpatrick 
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ground 


0.52 
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ground 
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ground 


2.16 


6.1 


1.7(0.4-2.8) 


0.1 (0.0-0.1) 


Host Fitzpatrick 


4% 


ground 


2.16 


6.1 
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ground 
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0.52 


8.0 
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Fig. 8. The fitted values of Rv with the Icr errorbars for all images that could be fitted using the Fitzpatrick parameterisation of dust 
in an intervening galaxy. The first panel shows Ry as a function of dust redshift and the second as a function of colour excess. In the 
first panel, the results from the different images of the same quasar will be plotted on top of each other since they all have the same 
redshift. In the second panel, the results of different images will be plotted in different places along the x axis since E{B - V) differs 
at least slightly between the images. The fitted values which have a larger probability than 10% (see TableO are plotted with stars 
and dark blue solid errorbars while the less reliable ones, with a smaller probability than 10%, are plotted with plus signs and dotted 
light blue eiTorbars. The most probable value of Ry found when co-adding the probability distributions for the fits (see Figure |7]l is 
indicated with a dotted line at Ry - 2.4. 



also b een investigated using the differential method by other ati - 
thors (iFalco et al.lll999HElfasd6ttir et al.ll2006tlToft et al.ll2000h . 
When we state their prefeiTed values of Ry in the text below fol- 
lowed by an interval within parenthesis, this is the error interval 
of Ry provided in that paper. 



SDSS J 131058.13-h010822 2 This quasar-galaxy system was 
found in lOstman et al.l ( l2006l) . by matching the coordinates of 
quasars with those of galaxies. The colours of the quasar indi- 
cate that it is most likely aff'ected by dust extinction. In that pa- 
per, it was found to have Ry - 3.4(1.7 - 5.7) using the CCM 
parameterisation, and Ry = 3.4(2.4 - 4.7) using the Fitzpatrick 
parameterisation. This agrees well with what was found in this 
analysis, where the preferred value of Ry was 3.5. 

Accordin g to the New Yo rk University Value-Added Galaxy 
Catalog (Bl anton et al.ll20()5l) . the foreground galaxy is a star- 
forming galaxy at z = 0.04. Fitting the quasar colours with 
starburst-like extinction according to the CAB law, we find a best 
fit E{B - V) of 0.3 (0.2-0.4), where the;^^^ provides a probability 
of 99%. 



Q2237-I-030 There are four images of this quasar, three of them 
can be fitted with dust extinction in the intervening galaxy, all 
with unrealistically low values of Ry and an E{B - V) between 
0.5 and 1.2, depending upon which image is studied and the pa- 
rameterisation that is used. The extremely low preferred values 
of Ry for image A and C, are difficult to justify physically due 
to the strange appearance of the extinction laws for low Ry (e.g. 
see the CCM parameterisation with Ry = 0.5 in Figure[Tll- Image 
D has a larger Ry value, with a higher chi-square probability. 
Image B, which has one deviating band, could not be fitted by 
dust extinction. 

The extinct i on pro perties of this system have been measured 
by iFalco et al] (1 19991) . using the differential method. They ob- 
tained a high value of Ry of 5.29(4.4 7 -6.11), which is inconsis- 
tent with our result. lEliasdottir et al.l (12006) used the same tech- 
nique and derived values of Ry between 2.9 and 3.1, with eiTors 
of the order of 1.5. In our analysis, values below 1.3, were ob- 
tained. 

We note th at the colours of the object have been shown to 
vary with time (iMoreau et al.ll2005l) and therefore there must be 
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an effect, in addition to or instead of dust extinction, affecting 
the colours, for example microlensing. This would explain why 
we derive best-fit values for Ry that are unrealistically low, and 
why completely different values have been obtained using other 
data sets. 



SDSS J1 14719.89+522923.1 This cas e of a quasar shin - 
ing through a galaxy, was discovered by Ostman et al. (2006). 
However, no Ry could be determined. Konig et al. (2006) have 
since detected absorption features from the intervening galaxy 
in the quasar spectrum strengthening the assumption that the 
colours of this quasar are affected by dust absorption in the in- 
tervening galaxy. In this paper, we found that the colours are 
best-fitted by dust in the intervening galaxy, with a low Ry. 

Acco rding to the New Yo rk University Value-Added Galaxy 
Catalog (iBlanton et al.l 12005 1). the foreground galaxy is a star- 
burst galaxy at z = 0.05. We find that Galactic-like dust describes 
the observed colours better than typical dust in starburst galax- 
ies. Using the CAB law for starburst extinction, we find a best-fit 
value of E{B - V) of 0.2, with ?iX^ probabiHty of 12%. 



SDSS J084957.97+51 082 9.0 This quasar with a foreground 
galaxy was discovered by Ostman et al. (2006). It was found to 
have a measured value of Ry of 1.7(0.7 - 2.9) using the CCM 
parameterisation, and of 2.2(1.5 - 2.9) using the Fitzpatrick pa- 
rameterisation. This agrees well with ourprefeiTed values of Ry 
for this system, of approximately 2.2-2.3. 

Acco rding to the New Yo rk University Value-Added Galaxy 
Catalog dBlanton et al J 120051) . the foreground galaxy is a star- 
burst galaxy at z = 0.07. We found that the observed colours 
of the quasar were reproduced more effectively using the 
Fitzpatrick extinction curve, than using an extinction curve that 
should be valid for starburst galaxies. However, the fit with 
starburst-like extinction according to the CAB law was also 
good. The best fit with E{B - V) - 0.5 provided a x^ proba- 
bility of 44%. 



SDSS1 1 55+6346 This quasar has two images. For one image, 
we were unable to explain the measured colours by dust ex- 
tinction, and for the other image, there is only a small proba- 
bility that the best-fit dust scenario is the true explanation of 
the colours. The magnitudes of image B could be fairly well- 
described by dust, if the observational errors were larger. Image 
A, on the other hand, has a deviating K band magnitude. This 
is one case where one could suspect problems with the ground- 
based observations. It would be interesting to observe this object 
further to investigate if the deviation is a result of the observa- 
tions, or if it is an interesting feature in the spectrum. 



Ry - 0.9(0.7 - 1.2), for the image we will refer to as image A. 
However, we did not manage to fit this image with intervening 
extinction here, and we have therefore no Ry with which to com- 
pare. Instead we found that it could be described by extinction in 
the host galaxy with an Ry of approximately 2.3-2.8. Image B, 
however, can be reproduced by assuming that there is dust either 
in the intervening galaxy or in the host galaxy. 



SDSS J0924+0219 There ai-e three images of SDSS 
J0924+0219. The colours within these images cannot be 
fitted by the extinction laws used in our analysis. 



CLASS B1 152 +1 9 9 Using the d i ffereri tial method, both 
iToft et alJ (I2OOOI) and lEliasd6ttir et all (l2006h have measured Ry. 
They obtained Ry - 1.3 - 2.1 and Ry - 2.1(2.0 - 2.2), respec- 
tively. In our analysis, we have not managed to find a good fit to 
the observed magnitudes for any of the two images. 



HE0435-1223 There are four images of this quasar. None of 
these images have measured colours that could be fitted using 
our extinction laws. 



Q0142-100 Using the differential method. iFalco et al.l (1 19991) 
and Eliasdottir et al. (2006) have measured Ry for this system. 
They obtained 3.11(2.11 - 4.11) and Ry = 4.7(4.0 - 5.4), re- 
spectively. We find that there is little extinction in this system, 
causing unreliable Ry determinations. 



HE0230-2130 None of the four images of HE0230-2 130 could 
be reproduced by taking into account the effects of dust, using 
both the HST, and the ground-based observations. However, if 
we use only the ground-based observations, some of the images 
can be explained by dust, in particular image A2. This is one 
of few images for which we measure a higher Ry, than for the 
Milky Way, 3.2 for the CCM parameterisation and 3.4 for the 
Fitzpatrick parameterisation. 



BRI0952-01 15 iFalco et al.1 d 19991) fitted a value of 3.10(2.10 - 
4. 10) for Ry, using the differential method. With intervening ex- 
tinction, we found lower values of Ry. 



WFI2033-4723 There are four images of this quasar. None of 
these images have colours that could be reproduced well by our 
extinction laws. 



MG1 654+1 346 The colours of MG1654+1346 could be ex- 
plained by SMC-like dust in the intervening galaxy, with a low 
Ry. However, the x^ of the fit is high. 

CXOCY J2201 32.8-3201 44 There are two images of CXOCY 
J220132. 8-320144, both of which could be explained by dust, in 
particular at the host galaxy with low values of Ry. 



SDSS J1 004+41 1 2 There are four images of this quasar. The 
observations could not be fit by assuming dust extinction for any 
of the images. The ground-based observations appear to follow 
a smooth curve, but the HST observations appear to be more 
scattered. However, using only the ground-based observations, 
we were unable to obtain a good fit, using the extinction laws 
used in this paper. 



SDSS J0903+5028 This object was studied bv lOstman et alj J1 004+1 229 There are two images of this quasar. We were un- 
( I2OO6I) . using a different data set. The CCM parameterisation able to reproduce the colours of either image, using our dust 
gave Ry - 0(0 - 0.7), and the Fitzpatrick parameterisation extinction laws. 
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MG041 4+0534 This is an interesting case of a quasar that is 
heavily reddened. There have been claims both that the red- 
dening is caused by dust in the lens galaxy, and in the host 
galaxy. iLawrence et alJ (Il995h reason that the most likely cause 
is dust extinction in the lens galaxy. They argue that (1) the 
observed spectrum can be well-explained by dust at interven- 
ing redshiftfl while extinction in the host galaxy is unable to 
reproduce the spectrum for Galactic extinction and barely for 
SMC-like extinction; (2) the separation between the light paths 
for the different images at the host is too small to explain the 
difference in colours between the images; in contrast, a sep- 
aration at the lensing galaxy would be able to reproduce the 
colour difference better; and (3) two of the reddest quasars that 
are known, MG 0414-H0534 and MG 1131+0456, both have 
an intervening galaxy, and it is unlikely that both have a fore- 
ground galaxy, and that the foreground galaxy is not respon- 
sible for the ext remely red colours. H owever, for the case of 
MG 11 3 1 +0456. iKochanek et"aLl jlobo.) concluded that dust is 
not responsible for the red colour, at either redshift. They pro- 
posed instead that stellar emission from the host galaxy, could 
explain the red colours. Tonry & Kochanek ( 1999) suggest that 
the red colours of MG0414+0534 are caused by dust extinction 
in the host. Their arguments are that ( 1 ) the visible arc of the 
host galaxy is bluer than the quasar images and that the extinc- 
tion must be uniform over a large region, because there is little 
difference in extinction between quasar images, and no notable 
difference in colour along the arc; and (2) that the colour and 
surface brightness of the lens galaxy agrees well with that of a 
passively-evolving early-type galaxy on the Fundamental Plane 
(lKeetonetaOl998i) . 

When fitting for dust extinction, we found that the colours 
of the quasar could be explained by presence of dust in the lens 
galaxy, but not by dust in the host galaxy alone. All images gave 
reasonably similar values of Ry, indicating that the foreground 
galaxy has a homogeneous dust content. Since there are strong 
claims of dust extinction in the host galaxy, we decided to mea- 
sure the likelihood of dust extinction in both galaxies. We fitted 
for dust extinction at both redshifts with one Ry for each red- 
shift, a value of E{B - V) for the host extinction, and a E{B - V) 
value for each line of sight through the foreground galaxy. We 
found that the best fit was given by no extinction in the host 
galaxy [E{B - V) = 0], and severe extinction in the foreground 
galaxy [Ry = 1.5, and E{B - V) = 1.6, 1.9, 1.5, 1.5 for the diff'er- 
ent sight lines through the galaxy]. Marginalising over E(B - V) 
for the host galaxy, we find that the x^ is almost as good for 
small extinctions E{B - Y) < 0.3 as for the best fit which had 
E(B - y) = 0. For higher values of E{B - V), however, the;^'^ is 
much worse (see the first panel of Figure|9]l. From our analysis, 
it appears that the bulk of the extinction occurs in the intervening 
galaxy, but we cannot exclude an additional small extinction in 
the host. 

The best fit for the intervening galaxy v&Ry ~ 1.5, regardless 
of whether we include extinction in the host or not. In the second 
panel of Figure |9] there is a marginalisation over Ry in the in- 
tervening galaxy, showing how \}s\q x^ changes for different Ry. 
Falco et alj(ll999|) found a similar value of /?,/. 1.47(1.32-1.62). 



Eliasddt tir et alj (|2006) found values ranging between 0.4 and 
3.8, depending on the images that they compared. 



SPSS J01 21 47.73+ 00271 8.7 The best-fit values for Ry from 
lOstmanetal] (|2006|) for this quasai- system (3.3 for the CCM 



parameterisation and 1.8 for the Fitzpatrick parameterisation), 
are within the 1 (j uncertainty of our fit. 

SPSS J1 459 07.1 9+002401 .2 This object was also studied by 
lOstman et alj (12006 ). In that paper, the CCM parameterisation 
provided aniJy value of 0(0-0.7), and the Fitzpatrick parameter- 
isation, an Ry value of 1 .8(0.4 - 3.5). These results are consistent 
with the values derived in this paper 



' They used the wrong lensing redshift, but according to our analysis 
with colour comparisons the claim is still valid for the correct redshift. 



SPSS J1 4461 2.98+0351 54.4 In lOstman et al.l (l2006h . the 
CCM parameterisation gave Ry - 0(0 - 0.7) and the Fitzpatrick 
parameterisation T^i/ = 1 •3( 1 .0 - 1 .6). In this paper, we found no 
dust model that could explain the colours in a satisfactory way. 
However, if we ignore the u band magnitude, the colours could 
be fitted by intervening dust with an Ry of 2.0. 

5. Summary and conclusions 

The reddening of quasars by foreground galaxies has been inves- 
tigated. At our disposal, we have had 21 quasar-galaxy systems 
with a total of 48 quasar images. By looking at the rest-frame 
B - y, we conclude that the foreground galaxy, close to the line 
of sight, must aff'ect the quasar colours. Using redshifted red- 
dened quasar templates, we fit for Ry and E{B - V) using the 
measured colours. 

Of the data that we consider, 22 images could be fitted 
by dust extinction. Assuming dust extinction in the intervening 
galaxy following the Fitzpatrick law, we find a most probable Ry 
value of 2.4. This value is lower than the Galactic mean value of 
3.1. However, the FWHM of our distribution of Ry was deter- 
mined to be 2.7 and thus the Galactic mean value is within the 
range of the spread. We note that our sample will not be repre- 
sentative of the universal galaxy population because the majority 
of our quasars were observed by strong lensing surveys, which 
are particularly efficient in finding massive galaxies. The large 
spread in values indicates that the dust content can vary signifi- 
cantly between different galaxies. However, we found no strong 
correlation between the values QiRy and redshift, colour excess 
or galaxy type. In the future, when a larger dataset with superior 
quality is available, a similar analysis could be of considerable 
benefit to, for example, the use of extinction-corrected super- 
novae as distance indicators. It would not be surprising if a cor- 
relation was found, because metallicity and elemental abundance 
ratios are expected to affect dust properties, and both quantities 
evolve with redshift and galaxy type. Furthermore, there are at 
least two major m echanis ms that create dust grains, the first is 
in AGB stars (e.g. Mathis 1990) and the second in core-collapse 
supernovae (e.g. Bank & Clavton 2003) . which could give rise 
to different types of dust. 

The uncertainties in the determination of the amount of dust 
are large, when taking into account the different possible dust 
models. This fact may be of importance for fl ux anomalies, used 
as a p robe of cold dark matter substructure (iKochanek & Dalai 
I2OOI . 

Another interesting observation is that when several images 
of the same quasar could be fitted by dust, the preferred values 
of Ry were similar. This could indicate that Ry has little spread 
within individual galaxies, but that it varies between galaxies. 
A caveat, however, is that for many of our quasars one or sev- 
eral of the images could not be fitted by the standard reddening 
laws we tried, potentially challenging the conclusion of homoge- 
neous dust properties within individual lensing galaxies, unless 
the colour outliers can be attributed to something different from 
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Fig. 9. The x^ values for different E{B - V) in the host galaxy (first panel) and Ry in the intervening galaxy (second panel) for 
MG0414+0534 . 



dust. In total, there were seven quasars in our sample that had 
several images, one or several of which could be fitted with dust. 
Three of those quasars had similar values of Ry, for all images, 
regardless of extinction law, two quasars had at least one image 
that could not be reproduced by assuming the presence of dust, 
one quasar had images where the Ry values did not agree with 
each other, and one quasar was consistent with experiencing no 
dust extinction. 

Several systems were found with low values of Ry, compa- 
rable to what has been measured from global fits, using Type la 
supernova data. However, the distribution of Ry for our fits sug- 
gests that a large fraction of the galaxies considered in this study 
are compatible with Milky Way dust. Due to the limited number 
of objects in our study, as well as the possibility of selection ef- 
fects and analysis bias, no strong conclusion can be drawn as to 
whether the extinction properties in foreground galaxies, along 
the line of sight to quasars, differ or not from extinction in host 
galaxies of Type la supernovae. 
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